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1 Introduction 


The Finite Element-Boundary Integral (FEl-BI) technique has been used to 
analyze the scattering and radiation properties of cavity-backed patch anten- 
nas recessed in a metallic groundplane. A program, CAVITY3D, was written 
and found to yield accurate results for large arrays without the usual high 
memory and computational demand associated with competing formulations. 
Recently, the FE-BI approach was extended to cavity-backed antennas re- 
cessed in an infinite, metallic circular cylinder. EXCALIBUR is a computer 
program written in the Radiation Laboratory of the University of Michigan 
which implements this formulation. 

This user manual will give a brief introduction to EXCALIBUR and some 
hints as to its proper use. As with all computational electromagnetics pro- 
grams (especially finite element programs), skilled use and best performance 
is only obtained through experience. However, we will comment on several 
important aspects of the program such as portability, geometry generation, 
interpretation of results and custom modification. 


2 Formulation 

EXCALIBUR implements the FEl-BI formulation for cavity-backed antennas 
recessed in metallic circular cylinders. The formulation imposes some re- 
strictions on the geometries which will be discussed. Principally, use of the 
BiConjugate Gradient-Fast Fourier Transform (BiCG-FFT) iterative solver 
requires uniform zoning on the aperture of the cavity. Thus, each surface 
patch has cylindrical-rectangular shape. The basis functions, dyadic Green’s 
function and associated field formulas are given in a paper concerning scat- 
tering by these structures [1]. Modifications of this formulation for antenna 
analysis was given in another paper [2]. 

The uniform zoning requirement causes some difficulty in modeling; how- 
ever, with some practice, these difficulties may be overcome. For example, 
the specification of the patch and cavity size must both be expressed by an 
integer number of edges (hence nodes). Thus, if the cavity is twice the size 
of the patch, one has no problem specifying the patch and the cavity with 
the same uniform grid. However, if the ratio of the patch and cavity sizes are 
not integers, discretization may not be possible. This is often the case with 
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a continuous wraparound cavity. Such a cavity is shown in figure 1 along 
with an example of a discrete wraparound array. If the cavity size and patch 
size are not convenient, you must either change the cavity size as possible or 
change the radius of the cylinder. If the radius is changed slightly, it will not 
effect the electromagnetic properties of the structure, but it may allow' uni- 
form discretization. Although the restrictions imposed by the uniform zoning 
requirement seems rather stringent, with practice, an antenna designer will 
find that EXCALIBUR is quite flexible. 

3 Compilation 

The first task in utilizing EXCALIBUR is to compile and link the various 
files of the program. The following files are required to run EXCALIBUR and 
are given on the distribution diskette: excalibuT.f, dyadic./^ matrixGenerate.f, 
preprocessor./, rcs.f, fft.f, gauss. inc, fft.inc and excalibur.inc. These files 
perform the following functions: 

• excalibur.f. Main program, BiCG solver, matrix building subroutines, 
FE-BI subroutine, impedance insert and various auxiliary subroutines. 

• dyadic./. Compute dyadic Green’s function terms for admittance ma- 
trix, on-surface and far-field Fock functions and gamma function. 

• matrixGenerate./. Boundary integral and FE matrix terms. 

• preprocessor./. Geometry/mesh generator. 

• res./. Radar Cross Section, far-zone dyadic Green’s function and plane 
wave excitation functions. 

• fft./. Forward and inverse Fast Tburier Transform subroutines. 

• gauss. inc: Numerical integration parameters for gaussian quadrature. 

• fft.inc: Include file for 2-D FFT subroutines. 

• excalibur.inc: Main memory allocation file also contains variable dic- 
tionary. 

The distribution disk also contains the following Fast Fourier Transform op- 
tions: 

• fftCooley./. Classic Cooley-Tukey Radix-2 algorithm. 

• fftSplitRadix./. Split-Radix algorithm by Sorenson [3]. 
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(a) 



(b) 

Figure 1: Illustration of two types of arrays: (a) wraparound array; (b) 
discrete cavity array 
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• fftCRAY.f: Calls the CRAY library vectorized, multitasked 2-D FFT 
routine. 

• fftCONVEX.f: Calls the veclib vectorized 2-D FFT routine. 

Prior to compilation, the user should copy one of these files into fft.f. The 
optimized routines for the CRAY and CONVEX architectures should be used 
when possible while in general the Split-Radix algorithm should give better 
performance than the Cooley-Tukey version. To enhance their efficiency, the 
Split-Radix and Cooley-Tukey files utilize a decimation-in-frequency forward 
transform and a decimation-in-time inverse transform in order to avoid the 
need to perform bit reversal. 

Another file which is included on the distribution diskette is convert- 
ToASCII.f. To save disk space, the geometry information is stored in a 
binary format by preprocessor. The program convertTo ASCII is included 
to produce a human-readable file (ASCII). The nodes, elements, edges, un- 
knowns and other useful information is provided in a easy to read (although 
disk space consuming) format. All the programs are compiled and linked by 
invoking the UNIX make utility. A Makefile has been provided on the distri- 
bution disk. To date, EXCALIBUR has been successfully compiled, linked 
and run on the following architectures/operating systems: SUN, DEC UNIX, 
HP 9000/Txx, IBM RS/6000, Silicon Graphics IRIS, CRAY and CONVEX. 
Three variable must be set within the Makefile: 

• FF: The Fortran compiler name for the architecture. 

• FOPT: The Fortran compiler options i.e. optimization, precision, etc. 

• LOPT: The name of any libraries required for linking. 

The user should uncomment these variables for the target architectures in 
Makefile. EXCALIBUR is constructed by simply typing make at the com- 
mand line, while the binary-to-ASCII conversion program is constructed by 
entering make convert. EXCALIBUR is invoked by typing excalibur at the 
command line while the conversion program is run with the command con- 
vertToASCII. Finally, the directory may be cleaned up of all object and 
executable file by typing make clean. 
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4 Geometry Information 

The binary geometry file created by preprocessor contains all the information 
concerning the physical structure under study except for the placement of any 
probe-feeds or lumped impedance posts. Therefore, it is important that the 
user be aware of the geometry entered into the FE-BI portion of the code. As 
previously mentions, the conrer/7bi45C7/program creates a human-readable 
file from the machine-readable geometry file. 

The first information provided in the resulting ASCII file is the header, 
which contains the number of nodes, number of edges, number of unknowns, 
etc. and an example of the header is shown in figure 2. The next field contains 
the node information. The information given is as follows (see figure 3): 




Column 1: Node number 


Column 2: 
Column 3: 
Column 3: 
Column 4: 
Column 5: 
Column 6: 


Radial (p) coordinate in centimeters. 

Angular {4>) coordinate in degrees. 

Axial (z) coordinate in centimeters. 

Layer number from top of the cylinder (aperture). 
Row number from lowest axial coordinate. 

Column number from smallest azimuthal coordinate. 


Each node is associated with a physical location (p,^,z) and a grid location 
(layer, row, column). 

Grid points must be used in the discretization of a geometry since the 
BiCG-FFT solver requires that each node pair lie an integer number of units 
apart. Thus, the distance between two nodes (primed and unprimed) on the 
surface of the cylinder is given by 

R{n,m]n\m) = — n')aA^+ {m — m')Az (1) 

EXCALIBUR distinguishes between grid points and nodes. A grid point can 
be thought of as the intersection of two lines of a piece of graph paper which 
is placed on the surface of the cylinder, A node is a grid point which lies 
within a cavity. The row and column number associated with a node actually 
is the row and column number of the grid point which formed the node. The 
first grid point which corresponds to the lower-left corner of the grid has row 
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Binary filename : test . exc 


ASCII filename: test .ascii 


NODE statistics: 

Total number of nodes: 50 

Number of nodes on the surface: 25 

Number of nodes along the metallic walls: 41 

Number of nodes on surface metallic patches: 0 

Number of nodes which are resistive: 0 


EDGE statistics: 

Total number of edges: 105 

Interior edges: 9 

Aperture edges: 24 

a) substrate edges: 24 

b) resistive edges: 0 

Metal edges (NOT unknowns) : 72 

I I I I I I I I I UNKNOWNS I I I I I I I I I — > 33 


ELEMENT statistics: 

Total number of elements: 16 

Surface Elements: 16 


Figure 2: Geometry header. 
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Node 

rho (cm) 

phi (deg) 

z (cm) 

layer 

row 

column 

1 

1.00000 

-5.00000 

-0.50000 

0 

0 

0 

2 

1.00000 

-2.50000 

-0.50000 

0 

0 

1 

3 

1.00000 

0.00000 

-0.50000 

0 

0 

2 

4 

1.00000 

2.50000 

-0.50000 

0 

0 

3 

5 

1.00000 

5.00000 

-0.50000 

0 

0 

4 

6 

1.00000 

-5.00000 

-0.25000 

0 

1 

0 

7 

1.00000 

-2.50000 

-0.25000 

0 

1 

1 

8 

i.ooooO 

0.00000 

-0.25000 

0 

1 

2 

9 

1.00000 

2.50000 

-0.25000 

0 

1 

3 

10 

1.00000 

5.00000 

-0.25000 

0 

1 

4 

11 

1.00000 

-5.00000 

0.00000 

0 

2 

0 

12 

1.00000 

-2.50000 

0.00000 

0 

2 

1 

13 

1 . 00000 

0.00000 

0.00000 

0 

2 

2 

14 

1.00000 

2.50000 

0.00000 

0 

2 

3 

15 

1.00000 

5.00000 

0.00000 

0 

2 

4 

16 

1.00000 

-5.00000 

0.25000 

0 

3 

0 

17 

1.00000 

-2.50000 

0.25000 

0 

3 

1 

18 

1.00000 

0.00000 

0.25000 

0 

3 

2 

19 

1.00000 

2.50000 

0.25000 

0 

3 

3 

20 

1.00000 

5.00000 

0.25000 

0 

3 

4 

21 

1.00000 

-5.00000 

0.50000 

0 

4 

0 

22 

1 . 00000 

-2.50000 

a.5DT)00 

0 

4 

1 

23 

1.00000 

0.00000 

0.50000 

0 

4 

2 

24 

1.00000 

2.50000 

0.50000 

0 

4 

3 

25 

1.00000 

5.00000 

0.50000 

0 

4 

4 

26 

0.90000 

-5.00000 

-0.50000 

1 

0 

0 

27 

0.90000 

-2.50000 

-0.50000 

1 

0 

1 

28 

0.90000 

0.00000 

-0.50000 

1 

0 

2 

29 

0.90000 

2.50000 

-0.50000 

1 

0 

3 

30 

0.90000 

5.00000 

-0.50000 

1 

0 

4 


<< Remainder of nodes truncated » 


Figure 3: Node information. 
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= 0 and column = 0. For a wraparound grid, the first grid point is physically 
located 4> = —180° and once again has row = 0 and column = 0. 

The next set of information provided is the edges which form the cavities. 
The information given is as follows (see figure 4) 


Column 1: 
Column 2; 
Column 3: 
Column 3: 
Column 4: 
Column 5: 


Edge number. 

Left (lower) node forming the edge. 

Right (upper) node forming edge. 

Unknown number (zero indicates a fixed edge (e.g. metal)). 
Orientation (^-, (f>- or z-directed). 

Type of edge (metal, substrate, resistive or interior). 


Each edge is associated with two nodes and hence has an orientation in the 
cylindrical coordinate system. If an edge is metal, since EXCALIBUR uses a 
total field formulation, that edge’s weight is fixed at zero. All other edges are 
unknowns which must be solved using the BiCG-FFT solver. A substrate or 
resistive edge is associated with the boundary integral while interior edges 
contribute only to the FE portion of the system. Currently, resistive cards 
are not implemented in EXCALIBUR. 

The next set of information related the unknowns on the aperture of the 
cavities to their edge number. It also includes the row and column number of 
that edge in the discretization grid. Although this information is useful for 
understanding the mechanics of the BiCG-FFT solver, it is of little interest 
to the general user. The given information is (see figure 5) 

• Column 1: Unknown number. 

• Column 2: Associated edge number. 

• Column 3: Row of this edge in the discretization. 

• Column 3: Column of this edge in the discretization. 


The edges which form each element of the mesh are given next. Each 
cylindrical shell element consists of eight nodes which form twelve edges. 
This information is useful in visualizing the mesh and could be hooked into 
a graphics package to generate a 3-D picture of the mesh. The prototype 
element is shown in figure 6 which displays the node numbering scheme. The 
information given by convertTo ASCII is (see figure 7) 

• Row 1, Column 1: Element number. 
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Edge 



Node 1 

Node 2 

Unknown 

Orientation 

Type 

1 

26 

1 

0 

rho-directed 

metal 

2 

27 

2 

0 

rho-directed 

metal 

3 

31 

6 

0 

rho-directed 

metal 

4 

32 

7 

25 

rho-directed 

interior 

5 

26 

27 

0 

phi-directed 

metal i 

6 

1 

2 

0 

phi-directed 

metal r 

7 

31 

32 

0 

phi-directed 

metal 

8 

6 

7 

1 

phi-directed 

substrate ' 

9 

26 

31 

0 

z-directed 

metal 

10 

1 

6 

0 

z-directed 

metal 

11 

27 

32 

0 

z-directed 

metal 

12 

2 

7 

13 

z-directed 

substrate 

13 

28 

3 

0 

rho-directed 

metal 

14 

33 

8 

26 

rho-directed 

interior L 

15 

27 

28 

0 

phi-directed 

metal 

16 

2 

3 

0 

phi-directed 

metal 

17 

32 

33 

0 

phi-directed 

metal 

18 

7 

8 

2 

ph i - d i r e c t e d 

substrate = 

19 

28 

33 

0 

z-directed 

metal 

20 

3 

8 

14 

z-directed 

substrate 

21 

29 

4 

0 

rho-directed 

metal 

22 

34 

9 

27 

rho-directed 

interior 

23 

28 

29 

0 

phi-directed 

metal 

24 

3 

4 

0 

phi-directed 

metal 

25 

33 

34 

0 

phi-directed 

metal 

26 

8 

9 

3 

phi-directed 

substrate 

27 

29 

34 

0 

z-directed 

metal 

28 

4 

9 

15 

z-directed 

substrate 

29 

30 

5 

0 

rho-directed 

metal 

30 

35 

10 

0 

rho-directed 

metal 


« Remaining edges truncated » 


Figure 4: Edge information. 
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Hill II III 


Unknown 

Edge 

Row 

Column 

1 

8 

2 

1 

2 

18 

2 

3 

3 

26 

2 

5 

4 

34 

2 

7 

5 

40 

4 

1 

6 

47 

4 

3 

7 

52 

4 

5 

8 

57 

4 

7 

9 

63 

6 

1 

10 

70 

6 

3 

11 

75 

6 

5 

12 

80 

6 

7 

13 

12 

1 

2 

14 

20 

1 

4 

15 

28 

1 

6 

16 

44 

3 

2 

17 

49 

3 

4 

18 

54 

3 

6 

19 

67 

5 

2 

20 

72 

5 

4 

21 

77 

5 

6 

22 

90 

7 

2 

23 

95 

7 

4 

24 

100 

7 

6 


Figure 5: Relationship between unknown number and edge number on aper- 
ture. 
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X 


Figure 6: Cylindrical shell element. 
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• Row 1, Column 2-5: /)-directed edges. 

• Row 2, Column 2-5: fli>-directed edges. 

• Row 3, Column 2-5: z-directed edges. 

The final set of information provided is the element parameters as shown 
in figure 6 which includes (see figure 8) 

• Row 1: Element number. 

• Row 2: Pa, Pb and t = pb — Pa- 

• Row 3: 4>i, (f>r and a = <f)r — (f>i- 

• Row 3: Zb, Zt and h = — Zb- 


5 Geometry Generation 

Having reviewed the geometry information provided by preprocessor through 
convertTo ASCII, we are prepared to generate some example geometries. 
Specifically, we shall look at radiation and scattering by a 2 cm x 3 cm 
patch antenna which is placed in a 5 cm x 6 cm cavity, a four element discrete 
array of such cavities and the same radiating array placed in a continuous 
wraparound cavity. These three variations of the same geometry exhibit the 
main classes of structures encountered in practice. Users will find it very 
helpful to check the entered geometry via convertToASCII prioT to running 
the solver part of EXCALIBUR. In particular, it is useful to discretize the 
cavity without patches present and retain the node information since it will 
be necessary to specify the row and column of the lower-left corner of each 
patch as well as the number of edges along each side. 

EXCALIBUR has a preprocessor module which generates the required 
mesh. It first generates the surface nodes which are then used to create 
the 3-D mesh. An in-house modeling package, such as SDRC IDEAS, may 
be used to create such a surface grid so long as the nodes are an integer 
number of units apart. It would be quite easy to interface such a package 
with EXCALIBUR by replacing the subroutine simp/eMesA with a universal 
file reader. However, we have found that the custom mesh routine provided 
with EXCALIBUR (simpleMesh) is sufficient for most modeling tasks while 
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Element Edges 


1 

1 

2 

3 

4 



5 

6 

7 

8 

_ 


9 

10 

11 

12 


2 

2 

13 

4 

14 



15 

16 

17 

18 



11 

12 

19 

20 

- 

3 

13 

21 

14 

22 

- 


23 

24 

25 

26 



19 

20 

27 

28 

r 

4 

21 

29 

22 

30 



31 

32 

33 

34 



27 

28 

35 

36 

= 

5 

3 

4 

37 

38 



7 

8 

39 

40 



41 

42 

43 

44 



<< Remaining elements truncated » 


Figure 7: Edges associated with each element. 
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Element size parameters: 

Element number: 1 


rhoA 

= 0.90000 

rhoB 

= 

1.00000 

t = 


0.10000 

phiL 

= -5.00000 

phiR 

= 

-2.50000 

alpha 

= 

2.50000 

zB 

= -0.50000 

zT 

= 

-0.25000 

h 

= 

0.25000 

Element 

rhoA 

number : 

= 0.90000 

rhoB 

2 

1.00000 

t = 


0.10000 

phiL 

= -2.50000 

phiR 

= 

0.00000 

alpha 

= 

2.50000 

zB 

= -0.50000 

zT 

= 

-0.25000 

h 

= 

0.25000 

Element 

rhoA 

number : 

= 0.90000 

rhoB 

3 

1.00000 

t = 


0.10000 

phiL 

= 0.00000 

phiR 

= 

2.50000 

alpha 


2.50000 

zB 

= -0.50000 

zT 

= 

-0.25000 

h 


0.25000 

Element 

rhoA 

number: 

= 0.90000 

rhoB 

4 

1.00000 

t = 


0.10000 

phiL 

= 2.50000 

phiR 

= 

5.00000 

alpha 

= 

2.50000 

zB 

= -0.50000 

zT 

= 

-0.25000 

h 

= 

0.25000 

Element 

rhoA 

number : 

= 0.90000 

rhoB 

5 

1.00000 

t = 


0.10000 

phiL 

= -5.00000 

phiR 

= 

-2.50000 

alpha 

= 

2.50000 

zB 

= -0.25000 

zT 

= 

0.00000 

h 


0.25000 


<< Remaining elements truncated » 


Figure 8: Element parameters. 
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being quite efficient. For this manual, we will use this package for all meshing 
requirements. 

The first example is a 2 cm x 3 cm patch antenna recessed in a 5 cm x 6 
cm cavity which is centered at (<;^ = 0®, z = 0 cm). The required information 
is as follows 

• Choose item 1 (Preprocessor) from main menu. 

• Enter radius of the cylinder in centimeters. 

• Enter angular and axial size of grid in degrees and centimeters. 

• Enter center of grid in degrees and centimeters. 

• Enter number of grid points in azimuthal and axial directions. 

• Enter number of identical cavities in the azimuthal and axial directions. 

• Enter number of nodes per cavity. 

• Indicate whether all surface nodes are metallic or resistive. 

• Indicate whether all surface nodes are on the substrate (0 means a 
patch is present). 

• Enter number of patches. 

• Enter row and column of lower-left hand node of the patch (see node 
section of the geometry file for this information). 

• Enter number of edges in (<?i,z) directions for this patch. 

• Indicate any additional metallic nodes (-999 2 denotes no remaining 
nodes). 

• Enter number of substrate layers. 

• For each layer, enter its thickness in centimeters. 

• Enter 1 to save this geometry. 

• Enter filename for this geometry. 

The transcript for this geometry is shown on the following pages. The 
next two examples are a discrete array of these elements and a continuous 
wraparound array. The main difference is that the number of grid points is 
not the same as the number of nodes on the cavity surface. Essentially, the 
same information is entered as above with the exception that these two arrays 
have a wraparound (360°) grid. See the attached transcripts for more de- 
tails. Note for continuous wraparound cavities, if a patch crosses the branch 
cut (^ = ±180°), the nodes along the lower edge of the patch must be 


76 


hand entered as shown in example 3. These nodes are obtained bj’ running 
the preprocessor without specifying any patches and inspecting the human- 
readable geometry file. Note that for these three examples, the number of 
unknowns is dramatically different. For the single cavity of example 1, only 
541 unknowns are required. The four identical discrete cavities requires 2164 
unknowns which is of course 4 x 541. Finally, the continuous wraparound 
cavity has 12992 unknowns! 
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I i I I I I I i M M i M I I M I I M II I I M ) I I i II I il I II t 

I I 

I Example 1 I 

{ Single Discrete Cavity Example I 

f I 

I I I II I I I I I I I I 1 I I I i I I I I I I I I I I I i I I I I I I I i I I I I 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 


Performing initialization, please wait.... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm); 

18.75 6 

Enter center of array (phi,z) in (deg, cm): 

0 0 

Enter number of grid points in (phi,z) direction; 
11 25 

Discretization : 

deltaPhi = 1.875 deg 

deltaZ » .25 cm 

Enter number of cavities in phi,z directions 
1 1 

Enter nxomber of nodes per cavity (phi,z): 

11 25 


275 surface nodes have been generated. . . 

The node numbering has the following pattern 

SINGLE CAVITY: 


1 

21 

22 

23 

24 

25 

1 

16 

17 

18 

19 

20 

z 

11 

12 

13 

14 

15 

1 

6 

7 

8 

9 

10 

1 

1 

2 

3 

4 

5 


if for example 25 nodes were specified, 
where node 13 is the center point. 


Now specify which other nodes are either: 
metallic 
or 

resistive 

All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l*yes) : 

0 

Enter number of metallic patches: 

1 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 1 
6 3 
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I 111 r 



Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

<< List of nodes constituting the patch deleted >> 

<< Each of these nodes is specified by (phi,z) coord >> 

Now enter any individual nodes... 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 

Also indicate if 2 = metal or 3 = resistive,.. 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered. . . 

66 metal nodes entered. . . 

Enter number of substrate layers (INTEGER) : 

1 


Generating sub-surface nodes... - 
Enter layer 1 thickness (cm) ; 
0.07874 

All 550 NODES have now been created 
which form 240 elements... 

Generating edges, .. 

1303 edges generated. . . , 


All element edges have been identified. . - 

Classifying edges. . . . 

Number of elements: 240 
Total number of nodes: 550 

Total number of edges: 1303 

Number of interior edges: 207 

Number of metal edges: 762 

Number of aperture edges: 334 

Resistive edges: 0 

Substrate edges: 334 

Number of unknowns: 541 

Determining unknown order. . . 

Do you want to save this geometry (1-yes) : 

1 

Enter EXCALIBUR filename: 
examplel . exc 

Writing binary EXCALIBUR file: 
examplel .exc 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3 ) Impedance aperture 
0) Exit 
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I M I I I i I M I M I M I II M i I I I I I I 1 I I I I M I i I i I 


j Example 2 

I Four Discrete Cavities Example 

I 

I I I 11 II I I I I M I I I I I I I I I I I I I I I I I I I I I I I I i I I I 
Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 


Performing initialization, please wait.... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm): 

360 6 

Wrap-around cavity entered. . . 

Enter center of array (z) in (cm) : 

0 

Enter number of grid points in (phi,z) direction: 
192 25 

Discretization: 

deltaPhi « 1.875 deg 

deltaZ * ,25 cm 

Enter 1 if ring element, 0 if rectangular array: 

0 

Enter number of cavities in phi direction 
4 

Enter number of nodes per cavity (phi,z): 

11 25 

Enter lower left-hand (row, column) : 1 
0 43 

Enter lower left-hand (row, column) : 2 
0 91 

Enter lower left-hand (row, column) : 3 
0 139 

Enter lower left-hand (row, column) : 4 
0 187 


1100 surface nodes have been generated. - . 

The node numbering has the following pattern 

AZMUTHAL CAVITY ARRAY: 


1 19 

20 

21 

22 

23 

24 

z 13 

14 

15 

16 

17 

18 

1 7 

e 

9 

10 

11 

12 

1 1 

2 

3 

phi - 

4 

5 

6 


if for example two 3x4 cavities 
were specified 

Now specify which other nodes are either: 
metallic 


or 
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resistive 


All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l=yes) : 

0 

Enter number of metallic patches: 

4 

Enter row, column of lower left corner of 
the patch on uniform grid for patch; 1 
6 46 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 2 
6 94 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 3 
6 142 

Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

« Node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 4 
€ 190 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Node list deleted » 

Now enter any individual nodes... 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 
Also indicate if 2 « metal or 3 - resistive.. - 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered, . . 

261 metal nodes entered. . . 

Enter number of substrate layers (INTEGER) ; 

1 


Generating sub-surface nodes... . 

Enter layer 1 thickness (cm) : 

0.07874 

All 2200 NODES have now been created 
which form 960 elements... 

Generating edges. . , 

5212 edges generated, . , . 


All element edges have been identified. . . 
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Classifying edges. . . , 



Number of elements: 960 
Total number of nodes: 2200 

Total number of edges: 5212 

Number of interior edges: 828 

Number of metal edges: 3048 

Number of aperture edges: 1336 

Resistive edges: 0 

Substrate edges: 1336 

Number of unknowns: 2164 

Determining unknown order, . . 

Do you want to save this geometry (l*yes) : 

1 

Enter EXCALIBUR filename; 
example2 , exc 

Writing binary EXCALIBUR file: 
example2 . exc 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 
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Example 3 

Four Patches on Continuous 
Wraparound Cavity Example 

I II I I I I I I I I I j I i I I i I I I I I I I I I I I I i I i I I I I I I I 


Do you wish to run: 

1 ) 

2 ) 

3) 

0 ) 

1 


Preprocessor 

FEM-BI 

Impedance aperture 
Exit 


Performing initialization, please wait.... 

Enter radius of cylinder (cm) : 

15.27887 

Enter array size (phi,z) in (deg, cm): 

360 6 

Wrap-around cavity entered. , . 

Enter center of array (z) in (cm) : 

0 

Enter number of grid points in (phi,z) direction 
192 25 

Discretization : 

deltaPhi * 1.875 deg 

deltaZ * .25 cm 

Enter 1 if ring element, 0 if rectangular array: 

1 

Enter number of identical rings (>=1) : 

1 Enter number of nodes per cavity (phi,z): 

192 25 

Enter lower left-hand (row, column) : 1 
0 0 


4800 surface nodes have been generated. . . 

The node numbering has the following pattern 

SINGLE CAVITY: 


1 

21 

22 

23 

24 

25 

1 

16 

17 

18 

19 

20 

z 

11 

12 

13 

14 

15 

1 

6 

7 

8 

9 

10 

1 

1 

2 

3 

4 

5 


if for exanple 25 nodes were specified, 
where node 13 is the center point. 


Now specify which other nodes are either: 
metallic 
or 

resistive 

All surface nodes metal or resistive (l=yes) : 

0 

Are the nodes substrate nodes (l~yes) : 

0 

Enter number of metallic patches: 
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4 


Enter row, column of lower left corner of 
the patch on uniform grid for patch: 1 
6 46 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

<< Metallic node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 2 
6 94 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

« Metallic node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 3 
6 142 

Enter number of unknowns (edges) for each 
dimension (phi, z) : 

4 12 

« Metallic node list deleted » 

Enter row, column of lower left corner of 
the patch on uniform grid for patch: 4 
6 190 

Enter number of unknowns (edges) for each 
dimension (phi,z): 

4 12 

This patch includes branch, , , 

Enter first row nodes: 

1343 

1343 176,25 -1.5 

1344 

1344 178,1249 -1.5 

1153 

1153 “180.0 “1,5 

1154 

1154 “178.125 -1.5 

1155 

« Metallic node list deleted » 

Now enter any individual nodes . . . 

Specify the class (metal or resistive) of all 
surface nodes which are NOT substrate nodes: 
Also indicate if 2 *= metal or 3 = resistive... 

Enter node number (-999 if done) : 

-999 2 

0 resitive nodes entered. . . 

261 metal nodes entered. , . 

Enter number of substrate layers (INTEGER) : 

1 


Generating sub-surface nodes.... 

Enter layer 1 thickness (cm) : 

0.07874 

All 9600 NODES have now been created 
which form 4584 elements... 

Generating edges... 

23616 edges generated. . . . 

All element edges have been identified... 


liil uiil I . , I i,|ii li i| 



Classifying edges.... 

Number of elements: 4608 
Total number of nodes: 9600 

Total number of edges: 23616 

Number of interior edges: 4416 

Number of metal edges: 10624 

Number of aperture edges: 8576 

Resistive edges: 0 

Substrate edges: 8576 

Number of unknowns: 12992 

Determining unknown order. . . 

Do you want to save this geometry (l^yes) : 

1 

Enter EXCALIBUR filename: 
examples . exc 

Writing binary EXCALIBUR file: 
examples . exc 

Do you wish to run: 

1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 
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6 Operation Modes 

EXCALIBUR has three main operation modes for FE-BI calculations (option 
2 from the main menu). They are: input impedance vs. frequency, radiation 
pattern or RCS vs. frequency and single frequency radiation and RCS pattern 
calculations. This section will describe each mode using example 1 above. 

6.1 Input Impedance - Multiple Frequencies 

The first option presented is calculation of a patch antenna’s input impedance 
at multiple frequencies. This is most useful in determining the resonant 
frequency of a patch antenna. The following information is required 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 1 to monitor convergence. 

• Enter 1 for diagonal preconditioning and 0 for no preconditioning. 

• Enter 1 for frequency sweep of the input impedance. 

• Enter name of file to store the input impedance. 

• Enter number of probe feeds. 

• Enter location of each feed in terms of ((^=degrees, z=cm). 

• Enter which layer in which the feed is embedded. 

• Enter complex current for this feed. 

• Enter number of impedance post loads. 

• Enter frequency range (in GHz) for this sweep. 

• Return to main menu. 

The input impedance for 3.1 GHz to 3.3 GHz computed every 10 MHz is 
shown in figure 9. The transcript for computing the input impedance of an 
axially polarized 2 cm x 3 cm patch in a 5 cm x 6 cm cavity is given on the 
following page. 
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Do you wish to run: 


1 ) Preprocessor 

2) FEM-BI 

3) Impedance aperture 
0) Exit 

2 

Performing initialization^ please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 

Material Parameter Specification. . . . 

Is the material filling constant (l=yeSr 0=no) ? 

1 

Enter relative permittivity [real, imaginary] : 

2.17 0 

Enter relative permeability [real, imaginary] : 

10 

Checking dimension allocations - 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 

Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (0=no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 = diagonal? 

1 

Do you want to compute: 

0 »= RCS/Pattern, 1 = Zin, 2 = Freq. Sweep? 

1 

Enter Zin filename: 
examplel . zin 

I I I I I II I I I I I I I I I I I I I I I II M I I I I I I I I I I I I I I 

I PROBE FEED INFORMATION 1 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

Enter number of feeds: 

1 

Enter feed location (phi,z) in (deg, cm): 

0 -0.375 

Enter layer number of feed <=; 1 

1 

Enter mag, phase of probe current (amp, deg) : 

1 0 

I I I I I I I II I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
I IMPEDANCE LOAD INFORMATION I 

I I I I I I I I I I I I I I I I I I II I I I I I I I I I I I I I I I I I I I I 

Enter number of loads: 

0 

Enter start, stop and increment freq. (GHz): 

3.1 3.3 0.01 

« Data shown in attached figure » 
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6.2 Pattern - Multiple Frequencies 

The next option presented is calculation of a radiation or RCS pattern at 
multiple frequencies. This is useful in computing the variation of gain or 
RCS with respect to frequency. Usually, a single observation angle is spec- 
ified although multiple angles are allowed. This mode permits radiation, 
bistatic and backscatter computations. In addition, the input impedance as 
a function of frequency is stored if a probe feed is used for excitation. In this 
example, we compute the backscatter at normal incidence for an £-;-polarized 
plane wave as a function of frequency. The required information is 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 1 to monitor convergence. 

• Enter 1 for diagonal preconditioning and 0 for no preconditioning. 

• Enter 2 for frequency sweep of the far-zone fields. 

• Enter name of file to store the input impedance. 

• Enter name of file to store the RCS or Gain. 

• Enter observation type (0 = backscatter). 

• Enter start, stop and increment azimuth (^) angles (in degrees). 

• Enter start, stop and increment elevation {B) angles (in degrees). 

• Enter polarization angle (0 = El-pol, 90 = H-pol). 

• Enter RCS filename. 

• Enter number of probe feeds. 

• Enter number of impedance post loads. 

• Enter frequency range (in GHz) for this sweep. 

• Return to main menu. 

The RCS for 3.1 GHz to 3.3 GHz computed every 10 MHz is shown in figure 
10. The transcript for backscatter calculations for backscatter calculations 
is shown on the next page. 
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Do you wish to run: 


1) Preprocessor , 

2) FEM-BI 

3) Impedance aperture 

0) Exit 

2 

Performing initialization, please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 

Material Parameter Specification, . , . 

Is the material filling constant (1-yes, 0=no) ? 

1 

Enter relative permittivity [real, imaginary ] : 

2.17 0 

Enter relative permeability [real, imaginary] : 

1 0 

Checking dimension allocations... 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 
Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (0=no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 * diagonal? 

1 

Do you want to compute: 

0 = RCS/Pattern, 1 = Zin, 2 = Freq. Sweep? 

2 

Enter Zin filename; 
test . zin 

Enter Frequency sweep data filename: 
examplel . f sw 
Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 


Enter start, stop, and increment azmuth angles [deg]: 

0 0 1 

Enter start , stop, and increment elevation angles [deg] 
90 90 1 

Enter polarization angle [0 <* alpha <= 90 deg] : 
Ephi=0 : alpha = 0 

Etheta=0: alpha = 90 

0 

Enter ouput RCS filename [<= 40 characters] : 
junk 

I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 
I PROBE FEED INFORMATION | 

I I I I I I I I I I 11 I I I I I I I I M I I I I I I I I I I I I I I I I I I 
Enter number of feeds: 

0 

I I I I I I I I I I I I I I I II I I I I I I I I I I I I I I I I I I I I I I I 
I IMPEDANCE LOAD INFORMATION | 

I I I I I I I I I I I I I I I I I I I I I M I I I I I I I I I I I II I I I I 
Enter number of loads: 

0 

Enter start, stop and increment freq. (GHz): 

3.1 3.3 0.01 


6.3 Pattern - Single Frequency 

The final operation mode is radiation and RCS pattern calculations at a 
single frequency. Of course, for the case of an antenna, the input impedance 
is also computed. This mode is generally used for multiple incident and 
observation angle applications. The only difference between this mode and 
the previous two modes is the option to save the admittance matrix in binary 
format. Since this matrix is excitation independent (though still frequency 
dependent), it can save time for multiple runs of a large geometry if this 
matrix is stored and then read in for each run. The required inputs are 

• Choose item 2 (FE-BI) from main menu. 

• Enter the stored binary geometry file. 

• Enter 1 if all elements have the same material parameters. 

• Enter complex permittivity. 

• Enter complex permeability. 

• Enter BiCG convergence tolerance, minimum and maximum number 
of iterations. 

• Enter 1 to monitor convergence. 

• Enter 1 for diagonal preconditioning and 0 for no preconditioning. 

• Enter 0 for single frequency operation. 

• Enter matrix storage/read option( 1 to store, 2 to read, 0 to do noth- 

ing)- 

• Enter six (6) character filename for matrix storage/read. 

• Enter observation type (1 = bistatic). 

• Enter incidence angle (<^,^). 

• Enter start, stop and increment azimuth {<f>) angles (in degrees). 

• Enter start, stop and increment elevation (0) angles (in degrees). 

• Enter polarization angle (0 = E-pol, 90 = H-pol). 

• Enter RCS filename. 

• Enter number of probe feeds. 

• Enter number of impedance post loads. 

• Enter frequency (in GHz). 

• Return to main menu. 

This option is illustrated in the next session transcript for bistatic scattering 
and radiation pattern calculations. 
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Do you wish to run: 


1) Preprocessor 

2) FEM-BI 

3) Impedance aperture 

0) Exit 

2 

Performing initialization, please wait.... 

Enter EXCALIBUR filename: 
examplel . exc 

Reading excalibur file: examplel.exe 
Material Parameter Specification. . . . 

Is the material filling constant (1-yes, 0=no) ? 

1 

Enter relative permittivity [real, imaginary] : 

2.17 0 

Enter relative permeability [real, imaginary) : 

1 0 

Checking dimension allocations... 

Number of Boundary Integral Unknowns: 334 
Total number of UNKNOWNS: 541 
Enter tolerance, minimum and maximum iterations: 

0.01 2 500 

Do you wish to monitor convergence (0=no, 1-yes) ? 

0 

Do you want: 0 = no preconditioning, 1 = diagonal? 

1 

Do you want to compute: 

0 = RCS/Pattern, 1 = Zin, 2 = Freq. Sweep? 

0 

Save the BI arrays in a file (0=*no, l=yes, 2=read) ? 

1 

All BI files will have the following form: 

<f ilename> . guu, <f ilename> . guv, etc . 

Enter BI binary filename (must be six (6) char): 
ex amp 1 

Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 

1 

Enter incident phi, theta [deg]: 

0 90 

Enter start, stop, and increment azmuth angles [deg]: 
-180 180 1 

Enter start , stop, and increment elevation angles [deg] 
90 90 1 

Enter polarization angle [0 <= alpha <= 90 deg] : 
Ephi=0: alpha = 0 

Etheta=0: alpha = 90 

0 

Enter ouput RCS filename [<= 40 characters] : 
exampl . bi 

Enter frequency (GHz) : 

3.3 

I I I I I I I I I I I I I I I I I I I I II I I I I I I I I I I I I I I I I I I 
I PROBE FEED INFORMATION I 

I II I I I I I I I M I I I I I I I I I I i I I I I I I I I I I I I I I I I I 

Enter number of feeds: 

1 

Enter feed location (phi,z) in (deg, cm): 

0 -0.375 

Enter layer number of feed <=: 1 

1 

Enter mag, phase of probe current (amp, deg): 

1 0 
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I I I i I I I I 1 { 1 I 1 M { I I I I I I I I I I I I I I 1! ! t I I I M I 
I IMPEDANCE LOAD INFORMATION I 

I I I I I I I I M I I I I II i I I I I I I I I I I I I I I I I I I i I I I I 

Enter number of loads: 

0 

Writing BI binary files... 

Generate bistatic excitation vector... 

Solve it . . , 

« Pattern deleted. » 

Run anouther excitation? (l=yes,0=no) 

1 

Enter observation type: 

0) Backscatter 

1) Bistatic 

2) Radiation 

2 

Enter start, stop^ and increment azmuth angles [deg]: 
-180 180 1 

Enter start , stop, and increment elevation angles [deg]: 
90 90 1 

Enter ouput RCS filename [<= 40 characters] : 
examp 1 . pat ■ 

Enter normalized pattern filename: 
exampl . norm 
Solve it . . - 

« Pattern deleted. » 

Run anouther excitation? (l^yes, 0=no) 

0 

Do you wish to run: 


1) 

Preprocessor 


2) 

FEM-BI 


3) 

Impedance aperture 


0) 

Exit 
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The bistatic patterns taken at 3.3 GHz for the three geometry configu- 
rations presented in this manual are shown in figure 11 The corresponding 
antenna pattern comparisonat 3.3 GHz is shown in figure 12 

7 Concluding Remarks 

This user manual presented some basic operation information for the FE- 
BI code, EXCALIBUR. This presentation was only meant to get an initial 
user started. As one becomes experienced with the code, additional features 
such as 2-D patch array modeling, multiple feed arrays and use of lumped 
impedance loads may prove useful. Indeed, an experienced user will find 
that custom features may readily be added to EXCALIBUR. For example, 
currently EXCALIBUR allows entry of material parameters either for the 
entire substrate, each layer of the substrate or on an element-by-element 
basis. This subroutine material \n file excalibur.f may readily be modified by 
the user to input a custom inhomogeneous substrate. 

The code is fairly “dummy proof”. If the user enters data which is not 
expected by EXCALIBUR such as a character when an integer is expected or 
an angle greater than 360°, the code will prompt the user to re-enter the re- 
quested data. Additionally, the storage allocation parameters in excalibur.inc 
must be set by the user prior to compilation. If a particular parameter is 
too small for a given run, the code will halt and suggest a new value for the 
offensive parameter. The user must reset that parameter, recompile and run 
the code again. The program also estimates the amount of RAM required 
at the start of a run. This estimate is based on the storage required by the 
arrays in excalibur.inc plus some scratch arrays. Each complex number is 
assumed to require eight bytes and each integer and real number require four 
bytes. The user should consider this estimate to be a slightly lower than the 
actual consumed memory. 
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Angle (<))) [deg] 


Figure 11: E- polarized bistatic patterns for a 2 cm x 3 cm embedded in the 
three different geometries presented in this manual. The incidence angle is 
(0°,90") and observation is in the 6 = 90® plane. 
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Angle (<(>) [deg] 


Figure 12: Axially polarized H-plane patterns for a 2 cm x 3 cm embedded 
in the three different geometries presented in this manual. 





Appendix 1 


A Note ox Applying New \’ector ABCs to Patch 
Antennas on a Circular Cylinder 


Leo C. Kempel and John L. Volakis 
Radiation Laboratory 

Department of Electrical Engineering and Computer Science 
University of Michigan, Ann Arbor, MI 48109-2122 

Recently, we developed a Finite Element-Boundary Integral (FE-BI) method which 
may be used to accurately model the scattering and radiation by conformal patch 
antennas. The use of a boundary integral to provide an exact mesh closure condition 
results in a rigorous formulation. When coupled with a Biconjugate Gradient-Fast 
Fourier Transform (BiCG-FFT) solver, the resulting computer code possesses both 
low memory and low computational demand. However, boundary integral closure 
involves a dyadic Green’s function which necessarily limits its extension to coated 
geometries. Although it is possible to construct a suitable Green’s function for mul- 
tilayered superstrates, the resulting formulation is both complex and bulky. Fur- 
thermore, this approach cannot be extended in a rigorous manner to more general 
doubly curved structures. 

Another hybrid FE method which has been used for very large 3-D scattering 
calculations at the University of Michigan is the finite element- absorbing boundary 
condition (FE-ABC) method. ABCs provide an approximate mesh closure which 
maintains a highly sparse matrix for the whole system. We propose utilizing a 
new conformal ABC (A. Chatterjee and J.L. Volakis, Microwave and Optical Tech. 
Letters^ 6, No. 16, pp. 886-889, Dec. 20 1993) which allows the closure surface to 
be brought quite close to the cavity aperture. As a result, the number of unknowns 
required is significantly less than the number required if a traditional spherical 
closure surface/boundary condition is used. 

We have applied the second order ABC to radiation and scattering by cavity- 
backed patch antennas which are embedded in a metallic cylinder. Although ex- 
tensive testing has not as of yet been accomplished, the preliminary results are 
encouraging. We looked at a 2 cm x 3 cm patch antenna which was printed on a 
5 cm X 6 cm X 0.07874 cm substrate which had a dielectric constant of = 2.17, 
The second order ABC was place 0.5A from the cavity aperture as shown in figure 
1. Two comparisons with the FE-BI formulation were performed. The first involved 
bistatic scattering with normal incidence {(f>i — 0®,^j = 90^) and observation in the 
0 = 90'’ plane for an E^-polarized plane wave which is shown in figure 2. A radiation 
pattern for the same antenna is shown in figure 3 where the resonant frequency is 
3.3 GHz. Since this is an axially polarized element, figure 3 is an H-plane pattern. 

The scattering and radiation example given herein involves a conformal antenna 
element without a covering superstrate. The FE-ABC formulation may be readily 
modified to account for such a covering which may consist of multiple inhomogeneous 
layers ets shown in figure 1. A future report will develop the formulation for both 
covered and uncovered elements and it will provide suitable validation of the method. 
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In addition, we shall investigate whether the near-field performance of the FE-ABC 
method is sufficient for input and mutual impedance calculations. Having shown 
that the FE-ABC approach can indeed be used for conformal antennas on singly 
curved surfaces, we may in confidence develop a similar formulation for antennas 
mounted on doubly curved surfaces. 



Figure 1: Coated cavity-backed patch antenna with ABC mesh termination. 
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Angle ((()) [deg] 

Figure 2: Comparison of 2"*^ order conformal ABC mesh closure condition 
with the exact boundary integral condition for bistatic scattering. The patch 
was 2 cm x 3 cm printed on a 5 cm x 6 cm x 0.07874 cm dielectric substrate 
(cr = 2.17). The ABC was placed 0.5A from the aperture and the operating 
frequency was 3.0 GHz. The £^^-polarized plane wave was incident from (<;i, = 
0®,^i = 90°) and observation is in the 6 = 90° plane. 



Angle (([)) [deg] 

Figure 3: Comparison of 2"'^ order conformal ABC mesh closure condition 
with the exact boundary integral condition for antenna pattern calculations. 
The geometry is identical to the one used in figure 2. The operating frequency 
wajs 3.3 GHz. The feed point was {(f)s = 0°,2, = —0.375 cm) which results in 
axial polarization and observation is made in the H-plane. 
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APPENDIX 2 


FE-BI implementation for coated planar platforms 

Consider the cavity recessed in a coated ground plane as shown in Fig 1. The 
coating is assumed to be of thickness d and of relative permittivity We 
will denote the free space region above the coating (z > d) as region 1, the 
region inside the coating (0 < z < d) as region 2 and that inside the cavity 
(— c < z < 0) as region 3. We will assume that the cavity is filled with 
an inhomogeneous material having a relative permittivity frc(r) and relative 
permeability /Xrc(r). In accordance with the equivalence principle, the fields 
in regions 2 and 3 can be decoupled by closing the aperture with a perfect 
conductor and introduce the equivalent magnetic current 

M = Exz (1) 

over the extent of the aperture, where E is the electric field at the aperture (z = 
0). The magnetic field is then due to the radiation caused by the equivalent 
current M residing on the ground plane and is given by 

H(r) = H‘(r) + H"(r)-jfcy[i+^Vv] -F (2) 

where F is the magnetic vector potential and is given in terms of the dyadic 
Green’s function of the coating as 

F = y^G(r,r')-M(r')d5" (3) 

and H* is the incident field, H’’ is that reflected by the coated ground plane 
without the aperture, S denotes the planar surface area of the aperture. If 
the conducting surfaces are restricted to horizontal planes, as is the situation 
here, we can find an equivalent scalar Green’s function Gv such that [1] 

V-G = V'Gv(r,r') (4) 


where 


g; agf _ /fcV Gf 

dz \ / J^y^T 


( 5 ) 
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More specifically, the function P defined in [2] is identically zero. Xoav consider 
Fs as 


Fs = -jk\ 

Making use of (4) in (6) we can write Fs as 

JY 


l+lvv 


F 


(G) 


Fs = -jkY JJ^G-M{r')dS' + ■M{r')Gv{T,r')dS’ (7) 

The required components of the Green’s function in (7) are given by 


G X3: — 


Gy.= 


Gyy = 


G.y = 


oo roo 




— OO J — oo 


X 


a oo 

-OO 




— OO J — oo 


'JP 1* A»y 


OO ^00 ^ 


— OO J —oo 


V^Mr A>y 


where 


Gxx = (er^o^yFi - k\klT-^ / {%ir^klkx^'^^ 

Gyy = -k\klV^) / 


( 8 ) 

( 9 ) 

( 10 ) 

( 11 ) 

(12) 

(13) 




-Kky (fc^Fa + CrklTi) / (STr^klkijS'^) 


with Tj and Fa given by 


Also, 


p _ e^fca sm{kid) — jki cos(feid) 
Cr^a cos(kid) + jki sin(A:id) 

p _ fea cos(fcid) + jki sin(fci</) 
k^ sin(fcid) — jki cos(kid) 

oi = r r 

J—oo •/— oo 


(14) 


(15) 

(16) 

(17) 
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with 


qh ^ r r ciu 

7-00 J-oo 


( 18 ) 


where 


and 


where 


[C cos{kiz) + D sin(tir)] 


ky [fci cos{kid) + jtrk 2 siu(fcid)] 
ki [cr^2 cos(fci(i) + jki sin(fc]d)] 



Gzx = 7^ [Esin(fciz) + Fcos(A:iz)] 


kx [ki sin(fcid) — jk 2 cos(A:id)] 
koZo [fci cos(fcid) + jk 2 sin(fcid)] 



(19) 

( 20 ) 
( 21 ) 

( 22 ) 

(23) 

(24) 


with k^ = erk^ - k\ = kl - and ^ kl + kl. After weighting (7) 
becomes 


Fs = -jkY M(r) • G • M(r')dS dS' 

IjT V . M(r) M{T')Gv{r, r')dS dS' (25) 

The fields in region 3 (the cavity region) are formulated employing the finite 
element method [3]. 


Evaluation of the Sommerfeld integrals using the Chebyshev de- 
composition method 

This is a new, fast and efficient method [4] of evaluating Sommerfeld integrals 
encountered, when we use a spectral-domain transformation to calculS'te fields 
of a dipole source in a homogeneous space or in a layered medium with planar 
boundaries. The integration is based on the decomposition of the integrand 
into Chebyshev polynomials, the oscillatory nature of which facilitates rapid 
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integration of the fast oscillating integrands. In (8)'(18) we encounter Fourier 
transfomi pairs given by 



1 f 

f{k., k.) = ^J dxdy 

(26) 


1 too t 

y) = ir / dk^ky 

Z7T J 7 — 00 

(27) 

Introducing the Bessel function as 



MKp) = ^ r 

Z7T Jo 



(28) 

we have 

» y*00 



f{kp) = / Jo{kpp)f{p)pdp 
Jo 

(29) 


foo 

fp= Jo{kpp)f{kp)kpdkp 

Jo 

(30) 


(29) and (30) is the Fourier- Bessel or Hankel integral transform pair. Inverse 
Hankel tranforms such as (30) are the generic form of Sommerfeld integrals 
which we consider for the demonstration of this niethod. Transverse derivatives 
of f{p) can be expressed as : 

d f 

- cos <l>J^ Ji{kpp)f{kp)kldkp (31) 

Hence the most general Sommerfeld integral is written as 

J fOO 

Jm{kpp)h{kp) dkp (32) 

0 

where is the Bessel function of the first kind of order m. Assuming that 
h{kp) has a pole at kp = kp {0 < kp < y/^ko) 

SI = / iJm{kpp)h(kp) - F{kp)) dkp 

Jo 

' F{kp) dkp 

0 

too 

+ / Jm{kpp)h{kp) dkp (33) 

J\/c^ko 
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with 


( 31 ) 


= 


f^p — kp 


where A is the residue of Jm{kpp)h{kp) at kp = kp. The second integral in (33) 
is solved analytically. Next we make the behaviour of the integrands in (33) 
smoother by applying a change of variables, for the first and third integrands 
respectively 

kp = ^/^kQ sin ( (35) 

and 

kp = y/^ko cosh ^ (36) 

Transforming the domain of integration to the domain of definition of the 
Chebyshev polynomials the first integral of (33) becomes 


If' 




(37) 


with 


/,,;( 2 ) = Re, Im {(J„(k,p)h{k,) - F{k,)),Jl - k^,VT^) (38) 

kp = v^fcosin (^{z + 1)) (39) 

Similarly transforming the domain of integration for the third integral in (33) 
and truncating the upper limit at a and b for the real and imaginary parts 
respectively we have 


2J-1 ,/TT 


0 


dz 2 -h j \ 


■ b f 4 ( 23 ) feizs) 




^3 




(40) 


with 


/a, 4 ( 2 ;) = Jm cosh 

(^y/^ko sinh ^ ~ 


( 41 ) 


fsfiiz) = Re, Im (^y/clkocosh 
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(43) 


Now f\{z) to fei^) an* fiited with Clirhvshov ])ol ynomials so that 


r ^ 


f{^) 


U=i 


Cl 

2 


with 




fc=i 


The orthogonal property of these polynomials is defined as 


/: 


' r,(z)r,(z) 


y/l — 


dz = 


0 (^7^ ;) 

f 

7T (z = j = 0) 


Applying the orthogonal property we evaluate the SI as 

^ 24(^14 +ic2,i) + 2 


with 


,2 - s(k)\ 

/(..c,)= g 


*+iCi,fc+i 


(44) 


(45) 


(46) 


(47) 


Using the Chebyshev decomposition method, a typical Sommerfeld integral 
encountered while analysing a horizontal electric dipole in the x direction 
above a ground with relative pemittivity and conductivity <r, is evaluated. 
The integral is expressed as 

foo /g-“oh+fc|\ 

A* = 2 / Jo{kpp) I 1 ) kpdkp (48) 

Jo \Uo + U J 

where 

P - UO = \jkl - P (49) 

u = Jkl - mV N^^er-j— (50) 

and k is the free space propagation constant. For this evaluation 2 (the ob- 
server height) to be equal to the dipole height h (A/50). The results of this 
evaluation is shown in Fig 2. 
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Region 1 



Region 2 Region 3 


Figure 1: Geometry of a groove in a coated plane 
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(V)sqv 


Variation of the Sommerfeld potential of a HED 

A 



®2 

Figure 2: A typical Sommerfeld integral evaluation 
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